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What is GMAT? 


■ GMAT is a mission design, 
optimization tool that is: 

- In-house 

- Open source 

- High fidelity 

■GMAT R2013a 

- Released April, 2013 

- 6 th public release 

1st major non-beta release 

■GMAT R2013b 

- Released August, 2013 

- Certification candidate 
Meets ACE requirements 
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What is GMAT?... 


■ GMAT can support flight regimes ranging from: 


- LEO 

- GEO 

- HEO 

- Libration 

- Lunar 

Interplanetary & Deep space 
- GMAT has supported 
- LCROSS 

- LRO 

- ARTEMIS 

- MAVEN 

- OSIRIS 

- TESS & more... 




Optimal Lunar Flyby 


Optimal Mars Trajectories 


Download and find out more: gmatcentral.org 



Outer Planet Transfers 



Asteroid (RQ36) Survey 
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ACE Mission Ov 


■ Sun-Earth LI Orbiter (Lissajous orbit) 

Spin stabilized & launched in August, 1997 

- Design amplitudes are: 

■ Ax = 81 ,755 km 

■ Ay = 264,071 km 

■ Az = 157,406 km 

- Sun-Earth-Vehicle (SEV) angle must be between 
4° & 20° nominal 


Lissajous orbit Viewed from Earth 



■ Station-keeping maneuvers: 

- 2-3 months apart 

- Nominal delta-V’s averaging 0.33 m/sec 

■ Attitude Maneuvers: 

- Performed weekly 

- Perturb ACE orbit 
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Flight Operation 



■ Began on August 2012 

■ Milestones 


- Requirements gathering 

- Gaps analysis for ACE requirements 

- Development/Testing/Documentation of new ACE related features 

- Develop ACE maneuver planning/product generation scripts & validate 
output 

- Write/perform Operational Procedures & Test Plans documents 

- Provide training to ACE Maneuver Team 

- Test Readiness Review 

- Non-Interfering Shadow Ops 

- Operational Readiness Review 
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Current ACE Od 


■ OD performed via GTDS 

■ Impulsive targeting/trajectory 
propagation performed via 
FreeFlyer 

■ Initial targeting done in ACE Eng. 
Coord, sys. Final targeting done 
in Attitude coord, sys. 

■ Finite-burn modeling is 
performed using GMAN 

■ GMAN generates Maneuver 
Cmd. File 

■ FreeFlyer delivers 28 days long 
ephem to NOAA 







Legend 

Tool 


Data File 

Manual Data Entry 






*Some information is manually entered and some is read directly from file 


We focused on tools/interfaces in red box 
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ACE Maneuve 


1. Get an updated OD state 

2. Prop to attitude re-orientation epoch & apply perturbations due to att. 
maneuver 

3. Next: Prop to maneuver epoch & enter Target Loop: 


Target Loop: 

Vary Z-component to 
Achieve RLP Vx = 0 @ 
RLP XZ plane crossing 
(i.e. When RLPY = 0) 
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Requirements G 


■ Requirements for ACE maneuver Ops gathered by working with 
maneuver planning team (97 requirements) 

■ Requirements had to be verifiable & unambiguous 

■ After 3.5 months of validation, final ACE Requirements approved 

■ ACE requirements areas: 

- Coordinate System 

- Force Model 

- Maneuver Targeting 

- Orbit Propagation 

- Product Output (SK dV, Code 500 & NOAA ephems & Maneuver summary report) 

- Spacecraft model 
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Gaps Analysis 


■ Missing features: 

1 . Parse through a vector hold file 

2. Write code-500 ephemeris file 

3. Develop new ACE Coordinate Sys. for maneuver targeting 

4. Report spacecraft acceleration 

■ 3 months of Development, Testing & Documentation efforts led to release 
of GMAT version R2013b (August, 2013)! 

■ R2013b is an internal release for Ops certification testing 
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Gaps Analyj 




Resources 


Output 


□■■Q Spacecraft 
DefaultSC 
h& Formations 
I ■■■■& Ground Station 
I ■■■■& Hardware 
0 -Q Burns 
frQ Propagators 
SolarSystem 
lil-C) Solvers 
i-fi) output 
E-3-& Interfaces 



h& Variables/Arrays/Strings 
0-Q Coordinate Systems 
!■■■■$. EarthMJZOOOEq 
!■■■■$. EarthMJZOOOEc 
EarthFixed 
L -^ EarthICRF 
H Functions 


© Mission | □ || B ' [■□■I 



T 

ii 

T 

=£ 

T 

Y 

& 




Filelnterface resource and Set command 


Code 500 ephemeris Format 


LocalAlignedConstrained Coord. Axis Type 


Spacecraft. ForceModel. Acceleration 
Spacecraft. ForceModel.AccelerationX 
Spacecraft. ForceModel. Acceleration Y 
Spacecraft. ForceModel. .Acceleration! 
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Pre-Shadow Ope 


■ There are two FreeFlyer scripts used for maneuver planning & product 
generation: 

- A CE_impulsive_ vec###. MissionPlan 

■ Generates weekly AV necessary to predict future SK maneuvers 

■ Used for both initial and final impulsive AV targeting 

- A CE_impulsive_NOAA 28day_ vec###. MissionPlan 

■ Generates 28 days long ephem. delivered to NOAA 

■ GMAT scripts were written using similar design philosophy: 

- ACE_impulsive_vec###.script 

- A CE_impulsive_NOAA 28day_ vec###. script 
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Local Operating Pr 



lopment 


■ ACE Maneuver team uses LOP document for End-to-End Ops support 
using FreeFlyer scripts 

■ Wrote detailed 45 page long LOP that instructs how to use GMAT scripts 
for ACE Ops: 

- Procedures for obtaining weekly ACE AV for Future Station-keeping Maneuver 

- Procedures for ACE Maneuver planning one week prior to the maneuver 

- Procedures for ACE Maneuver planning one day prior to maneuver 

- Procedures for final SK Maneuver planning (Post- Attitude Maneuver) 

- Procedures for generating NOAA 28-day Ephemeris 

- Procedures for delivering products via DataViewer 

■ Our LOP doc has been reviewed & approved by maneuver planning team 
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Test Plans Devel 


■ Wrote test plans for 97 requirements sub-divided in 6 areas: 

■ Coordinate System 

■ Force Model 

■ Maneuver Targeting 

■ Orbit Propagation 

■ Product Output 

■ Spacecraft model 

■ Each test plan: 

- Has detailed test procedures to test & verify each requirement 

- References separate GMAT ACE scripts to test each requirement 

■ACE team implemented test plans & GMAT passes all test plans & 
meets all requirements! 
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Requirements to 


Test Plans for Maneuver Targeting area: 


REQID Q 

Object Text B 

Test Plans 

MT01 

The ground system must be capable of ingesting the state vector from the TCOPS Vector Hold 
Files without user input. 

Follow procedure in FDSS-FORM-0015 Maneuver targeting Test 
Plan.docxand use ACE_impulsive_Burn_450.script GMAT script. 

MT02 

The ground system must be capable of ingesting the epoch from the TCOPS Vector Hold Files 
without user input. 

Follow procedure in FDSS-FORM-0015 Maneuver targeting Test 
Plan.docxand use ACE_impulsive_Burn_450.script GMAT script. 

MT03 

The ground system must be capable of ingesting C_r from the TCOPS Vector Hold Files 
without user input. 

Follow procedure in FDSS-FORM-0015 Maneuver targeting Test 
Plan.docxand use ACE_impulsive_Burn_450.script GMAT script. 

MT04 

The ground system shall use a user-input maneuver epoch for impulsive targeting. 

Follow procedure in FDSS-FORM-0015 Maneuver targeting Test 
Plan.docxand use ACE_impulsive_Burn_450.script GMAT script. 

MT05 

The ground system shall support varying the delta-V along the spacecraft body Z-axis during 
differential correction of impulsive maneuver targeting. 

Follow procedure in FDSS-FORM-0015 Maneuver targeting Test 
Plan.docxand use ACE_impulsive_Burn_450.script GMAT script. 

MT06 

Ground system shall propagate spacecraft to a user-specified number of XZ plane crossings 
in the Rotating Libration Point (RLP) frame during differential correction of impulsive 
maneuver targeting. 

Follow procedure in FDSS-FORM-0015 Maneuver targeting Test 
Plan.docxand use ACE_impulsive_Burn_450.script GMAT script. 

MT07 

The differential corrector shall compute a delta-V vector which achieves an accuracy better 
than 0.00000 ± 0.000001 km/s along X component of the velocity in the RLP frame (e.g., the 
Earth-Sun line) on the fourth X-Z plane crossing. 

Follow procedure in FDSS-FORM-0015 Maneuver targeting Test 
Plan.docxand use ACE_impulsive_Burn_450.script GMAT script. 
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Test Readiness 


■ On 09/10/2013: Presented TRR to ACE Ops Team 

■ Verify environment & tools are ready for shadow operations 

■ GMAT passed TRR! 


GMAT Development Team 

NASA GSFC 


16 



FDF Training 


■Although GMAT ACE LOP document serves as training & 
instructions manual to support ACE Ops, extra hands-on training 
was provided 

■ On 09/16/2013, gave training to maneuver team & demonstrated 
how to use: 

- GMAT’s ACE_impulsive_vec###.script 

- ACE_impulsive_NOAA28day_vec###.script 

■ Maneuver team now fully trained to use GMAT ACE maneuver 
planning & product generation scripts 
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Results/Analysis 

■ Delta-V comparisons 

■ Propagation comparisons 

■ Shadow Ops 

■ Operational Readiness Review 
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AV Comparisons 


■ SK AV validated against FreeFlyer using historical OD solutions 


AV diff. (perturbations from attitude re-orientation maneuver not modeled): 


TVHF file 

Maneuver Epoch 

GMAT SK AV 

AV diff. 


[UTCG] 

[cm/sec] 

[mm/sec] 

Vec424.txt 

15 Jan 2013 17:30:00.000 

15.01 

0.024 

Vec433.txt 

15 Apr 2013 16:00:00.000 

22.75 

0.017 

Vec440.txt 

19 Mar 2013 16:00:00.000 

12.53 

0.018 

Vec456.txt 

25 Jun 2013 19:15:00.000 

27.98 

0.016 


AV difference must be < 0.05 mm/sec 
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AV Comparisons 


AV diff. (perturbations from attitude maneuver modeled): 


TVHF file 

Att. Re-orientation 
Epoch [UTCG] 

Maneuver Epoch 
[UTCG] 

GMAT SK AV 
[cm/sec] 

AV diff. 
[mm/sec] 

Vec420 

19 Nov 2012 15:59:50.000 

19 Nov 201217:30:00.000 

29.65 

0.021 

Vec430 

15 Jan 2013 16:03:08.000 

15 Jan 201317:30:00.000 

19.97 

0.015 

Vec450 

02 Apr 201 3 17:49:36.899 

02 Apr 201 3 19:15:00.000 

19.47 

0.018 

Vec472 

09 Jul 2013 16:42:37.000 

09 Jul 2013 17:40:00.000 

15.30 

0.012 


AV difference must be < 0.05 mm/sec 
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Propagation Co 


■ Propagation compares using OD solutions from 4 TVHF files 

Short & Long term propagation comparison between GMAT & FreeFlyer: 


TVHF file used 

RSS position error after 28 
days in EarthMJ2000Eq 
[mm] 

RSS position error after 180 
days in EarthMJ2000Eq 
[meters] 

Vec433.txt 

0.50 

2.72 

Vec440.txt 

2.9 

3.04 

Vec450.txt 

6.1 

2.62 

Vec456.txt 

1.6 

4.73 

RSS pos. error (28 Days) must be < 10 mm 

RSS pos. error (180 Days) must be < 5 meters 
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Non-Interfering 


■On 09/23/2013, ACE maneuver team used GMAT & performed 
shadow operations during ACE SK maneuver: 

- Delivery products from GMAT verified against FreeFlyer 


AV diff. (perturbations from attitude re-orientation maneuver modeled): 


TVHF file 

Initial State Epoch 

AV diff. [mm/sec] 

RSS position error 

used 

[UTCG] 


after 28 days in 
EarthMJ2000Eq 




[mm] 


Vec493.txt 23 Sep 2013 00:00:00.000 0.015 1.83 

AV difference must be < 0.05 mm/s 
RSS pos. error (28 Days) must be < 10 mm 
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Operational Rea 


■On 11/19/2013: Presented ORRto ACE Maneuver Team 

■ Presented results from shadow Ops & test plans: 

- GMAT meets all requirements & passes all tests for ACE Maneuver 
Planning 

■ GMAT was deemed Flight Certified to support operational 
maneuver planning for ACE! 
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Conclusions/lm 



■ Demonstrated GMAT is flight quality software & is now Ops 
certified for ACE 

■ Laid groundwork for broad adoption of GMAT as an Ops tool for 
other GSFC missions 

■Goddard’s GMAT R2013b and recently R2014a: 

- Robust trajectory optimization tool available to all! 

- Provided a tool that Goddard controls to meet its unique and strategic 
needs 

- Provided a system for development of new mission design and nav. 
technology 

- In-house tool that complements other tools like FreeFlyer and STK 
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■ Requirements Gathering, 2001 

■ Architectural design, 2002 

■ Implementation of System Core, 2003 

■ First Beta Release, 2007 

■ Second Beta Release, 2008 

■ Decision to use as Primary Operational Software, 2010 

■ R201 la Release, 201 1 

■ R2012a Release, 2012 

■ R2013a Release, April 2013 (Production Release) 

■ R2013b Release, Aug 2013 (Ready for Ops Testing) 

■ Sep. 2013: NPR/GPR 7150.2 compliant 

■ R2014a Release, May 2014 
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Pre-Shadow Ops 


Basic Design methodology for GMAT’s ACE_impulsive_vec###.script : 


Resources | Mission | Outpu t | 


B t > Mis&on sequence 

■ user Inputs 
"tiff set I nits el Conditions 
B ♦ NO -burn Propagation fleg 

no Burn Propegetion - upto 2 weeks 
stop 1 

End No-burn Propagation fleg 
B ♦ if Apply Reorr dv 

j Apply Attitude re -orientation maneuver 
-«J End If Apply Reorr dv 
Prop to Maneuver Epoch 
Write prc-SK data 
B If spin^axis - 1 

iM Set Attitude from RA/DEC 
Else spin^axis — 1 
Set Attitude to ACE Coords 
■ 4 /M End If spin axra — 1 
US Loop Init 

R "Oj For I ^ lrnumCrossings 
® Init ACE S/C 
L-j. Target Loop 

m Vary MCC.V 
Set maneuver 
A V Apply MCC 

■ Clone s/c 

Decrement mass 
Bl iJOtl For J - IeI 

Prop to XZ Plane 
4 I End For J - 1st 
tef. Achieve VX — 0 
*J End Target Loop 
End For I — linumCrosangs 
Targeter Convergence Test 
Write post-SK. Data 
Hsr Compute total mass of ACE 

Prop to post Maneuver Epoch 
Create Ephem 


<=> 
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ACE Station-keeping 


■ 


r Context 


■ Initial Maneuver targeting is performed in ACE Engineering CS defined as 
follows: 

Z-axis: Defined by Earth center to ACE radial vector 

X-axis: up orthogonal to z-axis, in plane formed by z-axis & North Ecliptic Pole 
(NEP) 

Y-axis: Z cross X 


■ Final maneuver targeting is performed using Spin-axis Attitude CS once spin 
axis attitude is known prior to maneuver 

Z-axis: Defined by spin-axis attitude expressed in mean J2000 RA/DEC 
X-axis: Up orthogonal to z-axis, in plane formed by z-axis & NEP 
Y-axis: Z cross X 


■ Weekly spin-axis attitude re-orientation maneuvers perturb ACE orbit & 
perturbations modeled using Local Vertical Local Horizontal (LVLH) CS 
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GMAT ACE_ 

“User I 



t##. script 
Event 


© C:\Users\rqureshi\Desktop\ACE LOP Final Scipts\ACE_impulsive_450 - Copy.script - General Mission Analysis Tool (GMAT) 




File Edit Window Help 


[fc feBISl EJ ► II I 

Resources | Mission | Output] 


j Ha ! 1(5 | O W | Animation U O tiJ Lkl I GUI/Script Sync Status: | 


B--Q Mission Sequence 


y Set Initial Conditions 
B ♦ No-burn Propagation flag 

No Burn Propagation - up to 2 weeks 
# Stopl 

«J End No-bum Propagation flag 
If Apply Reorr dV 

y Apply Attitude re-orientation maneuver 
1 4 J End If Apply Reorr dV 
y Prop to Maneuver Epoch 
/ Write pre-SK data 
($•••♦ If spin_axis = 1 

Set Attitude from RA/DEC 
4* Else spin_axis = 1 
y Set Attitude to ACE Coords 
4 J End If spin_a»s = 1 

S y Looplnit 

S ® For I = l:numCrossings 
InitACES/C 
Target Loop 
!• Q) Vary MCC.V 
Set maneuver 
AV Apply MCC 

E Clone s/c 

Decrement mass 
For J = 1:1 

Prop to XZ Plane 
End For J = 1:1 
is4 Achieve VX = 0 
4 I End Target Loop 
J End For I = l:numCrossings 
y Targeter Convergence Test 
y Write post-SK Data 
y Compute total mass of ACE 
y Prop to post Maneuver Epoch 
y Create Ephem 


C) User Inputs 


0[f]0 


Initializations and user inputs 


% Input: state 

GMAT use_tvhf =1; % 1 = use tvhf, 0 = input state and Cr manually 

GMAT no_burn_prop = 0;% 1 = implement no-burn propagation, 0 = no-burn propagation turned off 
GMAT NumDays =14; % Enter # of days for no-burn propagation 

GMAT epochl = 'll Dec 2012 22:31:25.000'; 

GMAT X = 156666.73138825359; 

GMAT Y = -1412110.8264248; 

GMAT Z = -464550.53789868; 

GMAT Vx = 0.31253723305138; 

GMAT Vy = -0.043726652432292; 

GMAT Vz = 0.0074297160328779; 

GMAT Cr = 1.541788385; 

% Attitude Reorientation 

GMAT model_reorr =1; % 1 = model reorr, 0 = no reorr 

GMAT reorr_epoch = '02 Apr 2013 17:49:36.899'; %midpoint of reorr burn 
GMAT reorr_radial = -0.000005; %km/s 
GMAT reorr_tangential = 0.000006; %Jon/s 
GMAT reorr_normal = 0.000000; %km/s 

GMAT reorr_fuel_used = 0.067739; %kg (should be positive) 

% Maneuver 

GMAT maneuver_epoch = '02 Apr 2013 19:15:00.000'; 

GMAT maneuver_seed = 0.00019467808017750504; % initial dv guess, km/s 

GMAT target_crossing =4; % Nth crossing of XZ plane, where achieve condition is to be met 
GMAT tank_mass = 115.3644*0.45359237; % kg, total fuel remaining in all tanks 

GMAT sk_fuel_used = 0.0602; %kg (should be positive) 

GMAT define_spin_axis =1; % 1 = supply RA and Dec, 0 = build attitude based on radial vector 
GMAT Z_RA = 15.90; % deg 
GMAT Z_Dec = -6.50; % deg 

% Post maneuver epoch 

GMAT postmaneuver_epoch = '02 Apr 2013 20:00:00.000’; 

GMAT NumDays_after_maneuver = 180; % Enter # of days for which post maneuver ephemeris needs to be generated for. Variable is used in 'Create Ephem' ScriptEvent. 

%Remaining variable definitions 

GMAT drymass = ACE . DryMass ; 

GMAT ACE . FuelTankl . FuelMass = tank_mass; 

% Variable used for Targeter Convergence test 

GMAT targetTol = 0.000001; % Targeter shall compute dV vector with accuracy better than 0.000001 km/s along RLP Vx component on LAST X-Z plane crossing. 


□ H 


Zi 1 4pply 1 I Cancel [ 
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GMAT ACE G 


^ C:\Users\rqureshi\Desktop\ACE LOP Final Scipts\ACE_impulsive_450 - Copy.script - General Mission Analysis Tool (GMAT) 


L°iaj 


File Edit Window Help 


lafcliai<g| ia^ aii>- ii i 

j Res ou rces^ Mission [ Output 


BHCJ Mission Sequence 
• User Inputs 
yf Set Initial Conditions 
B ♦ No-bum Propagation flag 

iVf No Burn Propagation - upto 2 weeks 
# Stopl 

J End No-burn Propagation flag 
B ♦ If Apply ReorrdV 

i -j^f Apply Attitude re-orientation maneuver 
! 4 End If Apply ReorrdV 
Prop to Maneuver Epoch 
-tffi Write pre-SK data 
B ♦ If spin_axis = 1 

■jjjf Set Attitude from RA/DEC 
•••♦ Bse spin_axis = 1 
-jjf Set Attitude to ACE Coords 
<|J End If spin_axis = 1 
' Looplnit 

; For I = l:numCrossings 

a lnit ACE S/C 
Target Loop 
Qf Vary MCC.V 
1^ Set maneuver 
aV Apply MCC 
Clone s/c 
Decrement mass 
For J = 1:1 

Prop to XZ Plane 
End For J = 1:1 
Achieve VX = 0 
4 End Target Loop 
4 End For I = l:numCrossings 
■yf Targeter Convergence Test 
■g Write post-SK Data 
•Sj Compute total mass of ACE 
yt Pr °P to Post Maneuver Epoch 
y? Create Ephem 








^ »r fit „ * 1:08 PM| 

£ LMI 1 





1 ■" l 4>) * 9/29/20l| 
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GMAT ACE SK 




C:\Users\rqureshi\Desktop\ACE LOP Final Scipts\ACE_impulsive_450 - Copy.script - General Mission Analysis Tool (GMAT) 




File Edit Window Help 

hi' fr a i « i is* is I ► 

Resources Mission CMtiwt 


jj | 0 E§ | O j Animation IJ f~l U Lil GUI/Script Sync Status: 








iwH 

Q ||P^M 

£• II 

1 1 ^ 
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Pre-Shadow Ops..l 


Basic Design methodology for GMAT’s ACE_impulsive_NOAA28day_vec###.script : 




F1--Q Mission Sequence 
: User Inputs 

> * Set Initial Conditions 

Prop to Product_start epod~i 
h Write NOAA Ephem 


<=> 
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